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The accuracy of the WKB approximation when predicting the energy splitting of bound 
states in a double well potential is the main subject of this paper. The splitting of almost 
degenerate energy levels below the top of the barrier results from the tunneling and is thus 
supposed to be exponentially small. By using the standard WKB quantization we deduce an 
analytical formula for the energy splitting, which is the usual Landau formula with additional 
quantum corrections. We also examine the accuracy of our and Landau formula numerically 
for the case of the symmetric double well quartic potential. 

§1. Introduction 

Apart from the perturbative and variational methods the semiclassical or WKB 
approximation is the most widely used approximation in quantum mechanics to ob- 
tain the analytic expressions. The semiclassics serves not only as an approximative 
tool but it also helps us to picture and to understand the global behaviour of eigen- 
functions and energy spectra of Quantum systems which is particularly important 
in the context of quantum chaosS^'B. Since the systematic study of the accuracy 
of semiclassical approximation is a difficult task, it has been attempted for sirnple 
systems, where in a few cases even exact solutions may be worked outi^'&'il'&'Q'. 

In this paper we analyze the energy splitting of doublets in a generic one- 
dimensional double well potential. The splitting due to the tunneling is supposed to 
be exponentially small and is thus lost a-priori in the semiclassical approximation by 
expanding the quantum phase in the h, power series. By using the standard WKB 
technique we deduce an analytical formula for the energy splitting which is the usual 
Landau formula with additional quantum corrections and can be formally written 
as 

S 

AE = A exp 



h 



where S is the usual classical action inside the classically forbidden region (between 
the two turning points) and A is called the tunneling amplitude, which can be written 
as an h power series. This formula is based on a linear approximation of the potential 
near the turning points and so is inapplicable in the case of the square potentialsa^. 
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First we introduce the basic definitions, then we present the WKB spUtting formula 
based on the higher order semiclassical WKB expansion and finally we study its 
validity for the quartic potential. We shall demonstrate that our formula is indeed 
a significant improvement of the Landau approximation!). Also, we shall mention 
some potential applications in the field of molecular physics. 



§2. Basic formalism 



Let us consider a particle in a one-dimensional symmetric double-well potential. 
The stationary Schrodinger equation of the system reads as 

= ( - + Vix)y{x) = Ei^ix) . (2-1) 

The Sturm-Liouville theorem (see, for example 0^) ensures that for one-dimensional 
systems there are no degeneracies in the spectrum. Let ijji and "02 be two exact and 
almost degenerate eigenfunctions of the Schrodinger equation corresponding to the 
energy eigenvalues Ei and E2 lying below the top of the barrier. We write the 
eigenfunctions '01 and -02 in terms of the right and left localized functions 

Mx) = ^(Mx) + Mx)) and M-x) = ^(Mx) - Mx)) ■ (2-2) 

It is easy to show that Eq =< t/;Q{x)\H\ip(j{x) >=< ipQ{—x)\H\'ipo{—x) >= ^{Ei + 
E2) and 

2^2 

AE = E2-Ei = VoWo(O), (2-3) 

m 

(see&'il*). We shall remark that so far this is an almost exact starting formula to 
calculate the energv, splitting, since the error committed in the derivation is only 
exponentially smallQ^'tP. Please note that this quantity is always positive, because 
the tail of the right localized eigenfunction iPq{x) at x = has the same sign for 
-00(0) and its derivative V'o(O)- 



§3. Semiclassical method 



The WKB splitting formula is obtained by inserting the WKB approximant for 
the right localized function 0o which is determined by the WKB expansion of the 
Schrodinger equationEJ^. A generic eigenfunction ^ of the Schrodinger equation can 
always be written as 

V'(x) = exp Qcj(x)) , (3-1) 

where the phase (j{x) is a complex function that satisfies the Riccati differential 
equation 

'^'^{x) + {-)(j"{x) = 2m{E - V{x)) . (3-2) 
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The WKB expansion for the phase is given by 



oo *- 

<y{x) = Y.{-fa,{x). (3-3) 



fc=0 ^ 

Substituting (|3-3|) into ( |3-2D and comparing hke powers of fi gives the recursion 
relation (n > 0)™ 

n 

a'l = 2m{E-V{x)), ^ a'^aVfe + Vi = . (3-4) 

k=0 

The first five orders in the WKB expansion are given in the Appendix. 

In particular, if we call a and h the two turning points (let a < b) corresponding 
to the energy E, the right localized wavefunction ipQ for the forbidden region x < a 
that enters the splitting formula is given by 



C 

i^oix) = ^=exp 

v\p\ 



dx + deven + CTodd 



(3-5) 



where 



aeven = E W I ) ! ^2k{x) = a2fc(|p(x)|) = / ^(1^(01)^^ , (3-6) 



k=l 

with a2k{-\p\) = -£2ki\p\) 

and 



c^odd = ^^^"^^ ^2fc+i (|p(3;)|); 5-2fc+i(x) = o-2fc+i (|p(x)|) = / cr'2k+i {\p{0\)d^ 

k=l '' 

(3-7) 

with CT2fc+i(-|p|) = (T2fe+l( 



where p = \j2m{y{x) — E) and if we define ct^(x) = (T^(p(x)) = fk{p{x)) then 
'^fc(^) ~ ^A;(p(^)) = fkiPix)), for /fc being the same function of an appropriate 
argument in both cases and ak{\p{x)\) is considered as a family of functionals. In 
evaluating the integrals and we cannot integrate naively on the real axis, 
because such integrals are divergent, but must take a Dartiai derivative w.r.t. the 
energy of certain fundamental complex contour integralQ'. For potentials with only 
one minimum, thus having only two turning points, this has been proven by Robnik 
and RomanovskiE^^ 

The constant C is determined by the Kramers correspondence formulalll' and 
normalization condition for the right localized function 



which reads as 



\Mx)\^dx, (3-8) 



/' iexpf2^Vx = l, (3-9) 

Ja P V ?l / 
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that is approximately an integral of the semiclassical probability density | -00(3^)1^ 
over the classically allowed region a < x < bu> . 

In this way the splitting formula, up to the 5th order, after taking into account 
straightforward recursion relation (p-4p for or ct^ becomes 



AE 



h 

m 



b 1 



■ dx 



a P 



exp 



dx 



l + 2?i(5-2(0) + 27i^ I <5-f (0) + 



lafdx' 



laldx 



2 r ^dx 

'?4(0) + -ai(0)+ 2^2(0):^^^ 



2?i -CT^(O) + 2^2 (0)^4 (0) 



2^2^(0) 



V 



/g f dx 
rb 
J a 



Idx 



+ 2 



/g f 



(3-10) 



This formula is the usual Landau EP formula for the energy splitting (1st order in 
fi for the tunneling amplitude) with additional quantum corrections (up to the 5th 
order in fi for the tunneling arnplitude) . We note that higher-order WKB corrections 
quickly increase in complexity IP'i^ but, in principle, they can be calculated from the 
equation ( |3-4[) . It is important to stress that our splitting formula is good if the 
potential is sufficiently smooth so that the linear approximation and so the Kramers 
correspondence relations are valid near the turning points. 

The same splitting formula (3-10) can be derived using the semiclassical scatter- 
ing formalism and Kramer corresporidence rules, for example as expounded by Iyer 
and Willli^) and by Will and Guinnt3^. However, it is necessary to go back to the 
very first step in their formalism to get the result (3T0) of our problem. We have 
done this and confirmed, as mentioned above, that the result is the same. 



§4. Quartic potential 



As one of our test models we consider the symmetric double well quartic poten- 
tial, given by 

V{x) = -2Bx^ + Ax^ , (4-1) 

where the parameters A and B are related to the potential barrier Vq and to the 
position of the minimum xq by 



A ' 



xo 



By using the following reduced variables 



n 



eff 



Xo Vo ' - xo^/mVo ' 

the quantum Hamiltonian operator of the system can be written as 



H 



2 



2x2 ^ -4 ^ 



(4-2) 



(4-3) 



(4-4) 
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The exact spectrum of operator ([4-4| ) is calculated by numerical diagonalization 
of the Hamiltonian matrix, written in the basis of harmonic oscillator in QUADRU- 
PLE PRECISION number format (32 digits)!). 

The quantities <to(0), (T2(0), 174(0), | dx, ^ dx and ^^^^ dx that enters 
the semiclassical splitting formula can all be expressed in terms of complete elliptic 
integrals of the first and the second kind eP . 

Please note that all the integrands of these quantities are strongly divergent at 
the turning points, but all the expressions in (3T0) can be made convergent jand finjle 
by taking partial derivatives with respect to E of certain finite expressions 0) ' 13 'Q^' 13. 

In figure 1 we show the tunnel- 
ing amplitude ^ as a function of the 
mean energy E of all almost degen- 

"■^ ° 0.028 

erate pairs of quartic potential with 
hgff = 0.03. We compare the exact 0.026 
results with the semiclassical ones at 
1st (Landau), 3rd and 5th order in 
h. We can observe that, as expected, a 0.022 
the agreement is better at the energies 
for which the well known semiclassi- 
cal criterion is fulfilled. The criterion ^^^^ 
demands the absolute value of the de 
Broglie wavelength to be small com- ""le 
pared to the typical scale of the po- 
tential. This is not the case at ener- -1 -0.9 o.s 0.7 -os 0.5 0.4 0.3 0.2 0.1 
gies near the top and at the bottom of 

the potential barrier but the accuracy Fig- 1- Tunneling amplitude A vs mean energy E 
increases when going with the energy Quartic potential with h,ff = 0.03. 

away from those two values. In figure 

2 (3) we plot the tunneling amplitude j4 as a function fieff for the first (fourth) 
pair of almost degenerate consecutive energy levels. As shown also in table 1, the 
semiclassical results approach the exact ones by increasing the perturbative order in 
h. Note that at the 5th order in h the agreement with exact result is up to the 8th 
digit. 



t n,f = 0.0.3; Exact + 

t 

Landau A 

^ Semiclassical up to . * 

A 

A ^ Semiclassical up to h • X 



§5. Conclusions 



In this work we have taken the first step towards a systematic improvement of 
the Landau formulaB*, which is the semiclassical leading order energy level splitting 
formula for pairs of almost degenerate levels in double well potentials. We have 
developed the algorithm for the semiclassical fi expansion series to all orders for the 
tunneling amplitude A (of equation (|l-l|)), and thus calculated explicitly the quantum 
corrections up to the 5th order. We have compared the semiclassical predictions 
with the exact results obtained numerically, for the case of the quartic double well 
potential. Our approach is based on the usual WKB expansion in one-dimensional 
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potentials and so the calculation of higher corrections can in principle be continued 
by the same method, although the structure of higher terms increases in complexity 
very quickly. We have also checked what happens in cases where the assumption 
implicit in the Landau formula (namely the linearity of the potential around the 
turning points) is not satisfied as in the case of double square well potential: We 
get a different result even in the leading semiclassical order eP. We should stress that 
the Landau formula!^ is indeed quite good approximation since it always yields the 
correct order of magnitude (the exponential tunneling factor is always correct) and 
even the tunneling amplitude is correct within the 5-50 % . 

It is our goal to work out a more direct WKB approach to the solution of the 
multi-minima problem, by the contour integration technique, based on requirine the 
single valuedness of the eigenfunction, as has been done by Robnik and SalasnichB^'B 
in the case of a single well potential. This is our future project. 

Finally we should mention important applications e.g. in the domain of molec- 
ular physicsEP'Ej^'O). For example, the NH^ molecule can be described by a quasi 
one-dimensional potential V{x) as a function of the perpendicular distance x of the 
nitrogen atom from the plane, and as such it has the double well form, with 
the top of the potential barrier at x = 0. In order to calculate the energy level 
splitting of the doublets of vibrational modes one needs exactly our theory. An- 
other example is the torsional motion of C2H4 molecule, where again we encounter 
an effectively one-dimensional double well potential. In case of C2HQ molecule, we 
find three potential wells where the tunneling effects again determine the splittings 
of energy triplet levels and a generalization of our theory would give an improved 
estimate of the splittings. 
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Appendix 

The first 5 orders in the WKB expansion are: 
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p"'p' 13 99 



2 p5 ' 2 p6 



297 p^ 



) 



(5-1) 



32 



1 




111 




- 144 




+ 510 




) 



The derivatives a'^^-s can be obtained by replacing p{x) by p{x) in the expressions 
above. 
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Pair 1; Exact + 
Landau A 
Semiclassical up to * 
Semiclassical up to H^i X 



* 

A 



* 

A 



i 

A 



i 

A 



0.03 0.04 0.06 0.06 0.07 0.08 0.09 



Pair 4; Exact + 
Landau A 
Semiclassical up to h^^ * 
Semiclassical up to h^j X 



* 

A 



* 

A 



0.025^^^ ' ' ' ' ' ' 

0.03 0.04 0.06 0.06 0.07 0.08 0.09 0.1 
he, 



Fig. 2. Tunneling amplitude A vs heff for the first pair of almost degenerate consecutive energy 
levels (left). The same plot for the forth pair of almost degenerate consecutive energy levels 
(right). Quartic potential. 
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Table I. Energy splitting AE for all almost degenerate pairs at Tieff = 0.03. E^ is the exact mean 
energy, AEi^andau, ^Ej^a and AEf^s are the semiclassical results at 1st (Landau), 3rd and 5th 
order, respectively. In the last column of the upper table: exact energy splitting in units of 
mean level spacing. 



Pair number 


Es 


AEsxact 


AEn/{E7+^ - E?) 


1 


-0.9578013510838623 


9.232(2)E-28 


1.291(6)E-26 


2 


-0.8743363182686136 


6.533698E-25 


9.140541E-24 


3 


-0.7923227236203907 


2.197647599E-22 


3.074474586E-21 


4 


-0.7118393444104010 


4.66709587543E-20 


6.52919406381E-19 


5 


-0.6329771419061873 


7.0097983270055E-18 


9.8065981172742E-17 


6 


-0.5558426330893027 


7.901944841994295E-16 


1.10546971247045E-14 


7 


-0.4805627052945287 


6.920784581323967E-14 


9.68206928062435E-13 


8 


-0.4072917698084478 


4.806630946488390E-12 


6.72440144371474E-11 


9 


-0.3362229365299848 


2.675674688376366E-10 


3.74322699989559E-09 


10 


-0.2676066496781732 


1.196767019706533E-08 


1.67425832453064E-07 


11 


-0.2017846855780184 


4.273586943459416E-07 


5.97868122857108E-06 


12 


-0.1392610433553887 


1.195177462312803E-05 


1.67203455862220E-04 


13 


-8.088844930248588E-02 


2.500127371158503E-04 


3.49763905139656E-03 


14 


-2.855578436131249E-02 


3.352190557610379E-03 


4.68966210972998E-02 



AELandau{Es) 


AE^.{Es) 


AErAEs) 


8.56386072023299E- 


28 


9 


33424576483334E-28 


9 


16991849898935E- 


28 


6.34511775769322E- 


25 


6 


53796204696428E-25 


6 


53332080224002E- 


25 


2.15758453279967E- 


22 


2 


19798578660872E-22 


2 


19763766003663E- 


22 


4.60423712282166E- 


20 


4 


66737049838592E-20 


4 


66709238092202E- 


20 


6.93436277656352E- 


18 


7 


01000176474235E-18 


7 


00979708264343E- 


18 


7.83077639810467E- 


16 


7 


902O7840692513E-16 


7 


90194445833250E- 


16 


6.86703907509728E- 


14 


6 


92086233123417E-14 


6 


92078451430764E- 


14 


4.77381431395796E- 


12 


4 


80667193943203E-12 


4 


80663098878842E- 


12 


2.65943820818162E- 


10 


2 


67569528644699E-10 


2 


67567477543535E- 


10 


1.19031027734624E- 


08 


1 


19677788146677E-08 


1 


19676714995564E- 


08 


4.25353592752249E- 


07 


4 


27365612082292E-07 


4 


27358920983707E- 


07 


1.19079278916304E- 


05 


1 


19524083756918E-05 


1 


19518371126262E- 


05 


2.49835207740361E- 


04 


2 


50142635772383E-04 


2 


50058564705645E- 


04 


3.46270869995442E- 


03 


3 


39489990444765E-03 


3 


33278741702504E- 


03 



